Disorder and dephasing effect on electron transport through conjugated molecular 

wires in molecular junctions 
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Understanding electron transport processes in molecular wires connected between contacts is a 
central focus in the field of molecular electronics. Especially, the dephasing effect causing tunneling- 
to-hopping transition has great importance from both applicational and fundamental points of 
view. We analyzed coherent and incoherent electron transmission through conjugated molecular 
wires by means of density- functional tight-binding theory within the D'Amato-Pastawski model. 
Our approach can study explicitly the structure/transport relationship in molecular junctions in 
a dephasing environmental condition using only single dephasing parameter. We investigated the 
length dependence and the influence of thermal fluctuations on transport and reproduced the well- 
known tunneling-to-hopping transition. This approach will be a powerful tool for the interpretation 
of recent conductance measurements of molecular wires. 



INTRODUCTION 

The study of electron transport through nanostruc- 
tures has been a primary interest in the field of molec- 
ular electronics in the last few decades. Recently, in- 
creasing attention has been paid to the main role played 
by phase-breaking processes in molecular junctions. [l|- 
l5j The transition from coherent (tunneling) to incoher- 
ent (hopping) regimes is one example of such a process 
which provides helpful information for the fundamental 
understanding of transport mechanisms and for appli- 
cations such as organic semiconductors. Such transitions 
have commonly been reported in many studies of electron 
transfer in molecular systems and electron transport in 
molecular junctions, for instance, self-assembled mono- 
layers, [H H polymers, and macromolecules such 
as proteins, [12[ DNAs, [13[ and organic semiconductors. 
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For further development of molecular devices, the in- 
trinsic mechanisms of charge transport including the de- 
phasing effects and influential factors such as thermal 
fluctuations need to be understood in detail. In order 
to evaluate the electronic conductance of nanostructures, 
the Landauer formula is commonly used. However, this 
method is restricted to address systems within the coher- 
ent regime. In principle, this method must be extended 
before it can be applied to organic semiconductors or dis- 
ordered polymers where the dephasing effects play an es- 
sential role in electron transport. Normally, many molec- 
ular aggregates widely used in organic semiconductor de- 
vices exhibit high disorder features. A proper theoretical 
description of the transport properties of such defective 
structures has to take into account not only band-like 
coherence but also dephasing mechanisms via a hopping 
representation. [14 1 



One of the main sources of dephasing processes is 
the interaction between electronic states and molecu- 
lar vibrations, referred as electron-phonon (e-ph) cou- 
pling. The e-ph interactions can be experimentally in- 
ferred from inelastic tunneling spectroscopy or from ther- 
mal resonance broadening. The challenge, then, concen- 
trates in incorporating e-ph interactions into the Lan- 
dauer approach. The current state of the art as followed 
by several theoretical works involves the use of density 
functional theory (DFT) within nonequilibrium Green's 
function formalism where e-ph interaction is modeled by 
self-energy terms that can be solved self-consistently. [15| 
However, such a strategy can be computationally de- 
manding. Such limitations can compromise the under- 
standing of the main role played by dephasing effects 
on the charge transport properties of molecular devices, 
leading to the search for less expensive models that at- 
tend to more qualitative interpretations. One example 
of such a model is Biittiker's scattering approach, [l], [l6| 
in which e-ph coupling is replaced by phenomenological 
voltage probes capable of inducing phase-breaking pro- 
cess into the molecular conductor. 

Other approaches using reduced density matrix ele- 
ments where their time evolution can be described by 
Redfield equation have also been adopted to describe the 
physics of electron transfer/transport in molecular sys- 
tems and their crossover trends between coherent and 
incoherent regimes. fl7l - fl9j | Nevertheless, it is highly 
desired to investigate such a prominent physical phe- 
nomenon using a parameter-free (or using as few parame- 
ters as possible) first principles methodology since a more 
detailed quantitative study can be rendered using such 
techniques. Such tecniques would enable us to investi- 
gate the relationship between the conformational change 
of a realistic molecule and the change of electron trans- 
port. 
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In this letter, we analyzed the main role played 
by dephasing effects in the transport behaviour of or- 
ganic polymer systems using the Hamiltonian version of 
Biittiker's scattering approach which is the D'Amato- 
Pastawski model. [20j We show that the conductance 
of the molecules can be dramatically modified due to the 
dephasing effect depending on the size of the molecules. 
Important fundamental questions such as the tunneling- 
to-hopping transition and influence of thermal fluctua- 
tion on transport are also addressed in this work. We 
demonstrate that dephasing has a significant impact in 
driving the transport regime of molecular systems from 
a quantum mechanical tunneling to Ohmic behaviour for 
longer molecular wires. 
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FIG. 1. (a) A schematic model of a molecular junction 
connected between two electrodes. Biittiker probes are in- 
troduced to simulate dephasing events, (b) PPV- and (c) 
pyrrole-oligomer is coupled to electrodes consisted of graphene 
nanoribbons (GNRs) via peptide linkers. 



coupling between source/drain electrodes and the molec- 
ular wire. The electronic propagator for the coupled 
system is represented by a retarded Green's function 
defined as G R {E) = [(E + irf)I - H M - E L - £r] -1 , 
where El/ r are the self-energy elements which include 
the influence of the contacts. The conductance of a ID 
molecule at low bias in the coherent regime is then ob- 
tained via Landauer's formula, G = G Tr[G R T L G A r R ] 
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where Go is the quantum conductance unit Go = 2^- 
and Tl/r represents the broadening function given by 
r L/R (£) = i[£ L/R (£) - 4 /R (£)]. 

To include dephasing events in the model, the 
D'Amato-Pastawski formalism is adopted. Within this 
picture, one includes additional self-energy terms that 
account for the dephasing processes. This function is 
defined as Gf s (E) = [(E + i V )I - H M - S L - £ R - 
EIHi^b.ii]" 1 , where N is the number of reservoirs 
and Ee.n are the self-energies due to Biittiker probes; 
Eb,,i = z7 e _ p /i/2 that is written in term of a dephasing 
parameter, j e -ph, associated to the e-ph coupling. 

The effective conductance is given by T e g = Tl. r + 
El,„ = i T Ltm W~l n T n>R , where T L3 describes the coher- 
ent contribution and the subsequent term takes into ac- 
count the incoherent events. IP is a Markov matrix 
formed by transmission T mn between the m-th and n- 



th probes, and the reflection i?,- 
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broadening caused by dephasing processes is defined as 
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which is related to the 



,m ^B.m-I ~ Je-ph 

dephasing rate, k = T m /h. Although the dephasing 
strength "f e -ph can be evaluated by the Fermi Golden 
rule quantitatively from first principles, which are com- 
parable to ksT, [3, [22 1 it is also possible to examine the 
influence of the dephasing events on charge transport by 
changing the rate of dephasing events. 



RESULT AND DISCUSSION 



THEORETIACAL FRAMEWORK 

We calculated the conductance of molecular wires us- 
ing the Landauer formula written in terms of Green's 
functions. [2l[ Figure QTa) shows a schematic picture 
of the studied system which is composed of a two ter- 
minal conducting channel connected to a central scat- 
tering region. A sequence of local dephasing probes is 
incorporated to describe incoherent events that affect 
the transport properties of the system. The electronic 
structure of the molecular wires is represented by a tight- 
binding Hamiltonian; H = + Vl + Hm + Vr + ff R , 
where i?L/R and Hm represent the left/right electrode 
and central contributions, respectively. Vl/r defines the 



Firstly, it is important to obtain the electronic trans- 
mission within the coherent regime in order to compare 
it to the incoherent cases. We chose two sample sys- 
tems for such analysis: PPV- [Fig. [ljb)] and pyrrole- 
oligomer [Fig. HJc)] polymers. At first, we use a simple 
tigh-binding model, which only considers the nearest- 
neighboring 7r-orbitals in order to capture the trend of 
the transmission profiles. The coupling strength of C-C 
bonds on the GNR electrodes [23| are simply parameter- 
ized as (3 = 2.66 eV. In the scattering region, the transfer 
integrals for double and single bonds are set as /?d = 1-2/? 
and f3 s = 0.8/3, respectively, and the on-site energies for 
each atom are e 2p = 0.0/3 (C), 0.902,3 (N), and 1.902/3 



(O). [2J| The Fermi energy of the systems are set to 0.0/3 
for simplicity. The coherent transmission results are dis- 
played in Fig. [5]for (a) PPV- and (c) pyrrole-based molec- 
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FIG. 2. Electronic transmission calculated for PPV- (upper panels) and pyrrole- (lower panels) based molecular wires as a 
function of the molecular length within tight binding approximation. Results are displayed within coherent [(a) and (c)] and 
incoherent [(b) and (c)] regimes assuming also different dephasing parameters. 



ular wires considering different lengths. The transport re- 
sponse of the systems is characterized by a transmission 
valley centered around the Fermi energy (Ep — 0.0/3). 
Such valleys become more pronounced as the molecular 
length increases and this results in the well-known expo- 
nential decay of the conductance at the Fermi level as 
shown in the logarithm plots of Fig. [2] (b) and (d) (red 
lines). 

Dephasing effects are subsequently included in the 
model and we investigate how they affect the conduc- 
tance behaviour as a function of molecular length. Bear- 
ing in mind that the e -ph coupling is much weaker than 
the transfer integrals, 25| we finely changed the dephas- 
ing parameters between 0.5 and 20 m/3 and studied their 
influence on the charge transport as shown in Figs. [2](b) 
and (d). One can see that the transmission is consider- 
ably enhanced in comparison to the coherent tunneling 
as the length of the molecule increases. A clear deviation 
from exponential behaviour can be observed which be- 
comes more evident with increasing dephasing strength. 

More sophisticated methods such as the density func- 
tional based tight-binding (DFTB) approach [2 



27] 



were also employed to verify the robustness of such 
crossover behavior. We relaxed the molecular junctions 
and calculated the transmission in the same way as ref 
(li} . At first, we relaxed the isolated molecular wires and 
ID GNRs individually. Then we coupled the relaxed iso- 
lated molecular wires between the single unit cells of the 



relaxed ID GNR electrodes. We relaxed them using peri- 
odic boundary conditions fixing the atoms at the bound- 
aries. After the relaxation, we placed the relaxed cou- 
pled systems between the relaxed ID semi-infinite GNR 
electrodes and calculated the transmission probabilities 
using the same approach. 

The Fermi energy for the system was determined from 
the band calculation of the ID armchair GNRs which 
was determined to be -4.41 eV. Because of the steric re- 
pulsion between hydrogen atoms, the carbon-carbon dis- 
tances connecting the unit cells in the horizontal direc- 
tion is longer than that of others. Therefore the armchair 
of GNRs is slightly dimerized giving rise to the band 
gaps p9| around 5 eV, this can be seen as abrupt drop 
in the transmission functions in Fig. [3ja) and (c) . The 
transmission profiles and length dependence of conduc- 
tance at the Fermi energy with/ without dephasing effects 
calculated using DFTB method are presented in Fig. [3] 
A similar behavior as the one obtained via the simple 
tight binding approximation in Fig. [2] is observed. From 
the linear trend of coherent tunneling, we determined 
the damping constants for PPV, 0.24lA _1 , and pyrrole, 
0.308A -1 . The exponential decay associated with coher- 
ent transport is quenched as the dephasing effects get 
stronger. 

Once more, we verified that dephasing dictates the con- 
ducting response for longer ID systems. The change on 
the transport behavior from direct coherent tunneling to 
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FIG. 3. Electronic transmission calculated for PPV- (upper panels) and pyrrole- (lower panels) based molecular wires as a 
function of the molecular length within the DFTB approach. Results are displayed within coherent [(a) and (c)] and incoherent 
[(b) and (c)] regimes assuming also different dephasing parameters. 



a hopping mechanism can be evidently seen in Fig. |4] 
which shows the resistance as a function of length for 
both molecular structures calculated with the dephasing 
strength being 20 mcV. The system undergoes a crossover 
from the exponential to the Ohmic regime with the es- 
timated resistance per unit length for PPV and pyr- 
role being 283.9Mf2/nm and 39.3M£7/nm, respectively. 
Figure [5] shows the resistance per unit length in the 
Ohmic regime as a function of dephasing strength for 
both molecular wires. We can see that the resistance is 
decreased because of the assistance from the incoherent 
transport. 




Molecular length ~ > (nm) Molecular length — > (nm) 



FIG. 4. Resistance of (a) PPV- and (b) pyrrole-based 
oligomers as a function of molecular lengths. The fitting 
curves for coherent tunneling and hopping transport are 
shown in green and blue, respectively. Dephasing strength, 
"fe-pht was set to 20 meV 

Since we want to assess the influence of conformational 
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FIG. 5. Resistance of (a) PPV- and (b) pyrrole-based 
oligomers per unit length in the Ohmic regime as a function 
of dephasing strength. 

fluctuations on transport in realistic molecular wires un- 
der dephasing conditions, we performed molecular dy- 
namics (MD) simulations at room temperature using the 
DFTB method and calculated the effective conductance 
along the MD trajectories. Strictly speaking, the MD tra- 
jectory should be calculated under non-equilibrium con- 
dition and the velocity for the MD should include the 
effect of the heat produced tunneling current, however, 
we took the at equilibrium condition ignoring the heat 
produced by tunneling current for simplicity. At least 
from this calculation it is possible to examine the change 
of conductance due to the change of molecular structures 
crude approximation. 
The motion of the atoms during the simulations has 
been integrated using the standard Velocity Verlet algo- 
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rithm with a time step of 1.0 fs and a total time dura- 
tion of 10 ps. The influence of structural fluctuations on 
the conductance is investigated considering that only the 
atoms in the scattering region are allowed to move. After 
the MD simulation, 100 snapshots were extracted from 
the MD trajectories. Then we calculated the effective 
conductance for the 100 stational geometries in order to 
see the change of conductance due to the change of the 
coordinate. The e-ph coupling was set to the same as 
room temperature, 25.8 meV. |3| 
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FIG. 6. Electronic transmission at the Fermi energy within 
coherent and incoherent regimes as a function of length on 
the molecular dynamic flies for PPV-based (upper panel) and 
pyrrole-based molecular wires (lower panel). In the incoherent 
regime, dephasing strength was set to 25.8 meV. 

Figure |6] shows how the coherent and effective conduc- 
tances disperse along the MD trajectory with the molec- 
ular lengths. The dominance of hopping mechanisms for 
longer molecules is still preserved along the entire MD 
trajectory. The dephased transmission clearly deviates 
from the exponential trend of coherent tunneling. 

Another important feature between these two trans- 
port regimes is the degree of dispersion observed on the 
logT v.s. length plots in Fig. [5] Strong fluctuations are 
seen for the length dependence in the coherent tunneling 
regime whereas narrower distributions appear when de- 
phasing effects are considered. This issue can be clarified 
by analyzing the variation of the root-mean standard de- 
viation (RMSD) of molecular coordinates with respect to 
a reference coordinate, in this case, the initial one. 

Figure [7] presents the molecular length dependence of 
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FIG. 7. (a) Fluctuation (RMSD) of two types of molecular 
wires along MD pathways as a function of molecular lengths, 
(b) Standard deviation of logT for two types of molecular 
wires with/without dephasing effects. 



RMSD (left panel) as well as the standard deviation 
(right panel) of the transmission logarithm (logT) for 
the coherent and dephased regimes. RMSD results in 
Fig-BJa) point out that longer molecules potentially have 
more configurational degrees of freedom and, for this rea- 
son, RMSD increases with the size of the junction. Both 
types of molecules present approximately the same trends 
in the RMSD curves. The standard deviation of logT in 
Fig. [ZKb) highlights how strongly the coherent /dephased 
transport fluctuates along MD trajectory. We can see 
that coherent tunneling is more sensitive to the molec- 
ular disorder induced by the thermal fluctuations than 
dephased transport. Thermal fluctuations can disrupt 
the 7r-orbital conjugation of the molecules, localizing the 
states along the molecular wires. Such breaking of the 
delocalized 7r-orbitals, which bridge the electron tunnel- 
ing between two contacts, results in a large dispersion of 
coherent transmission in the logT v.s. length plots in 
Fig. E 

On the other hand, in the dephased regime, unlike in 
the coherent regime, dephased transport does not fluctu- 
ate strongly with the breaking of 7r-orbital derealization 
since there exists many individual jumps with probability 
T m ,n for sequential hopping processes and there always 
will be pathways for hopping that favor transport (note 
that dephased transport is determined by the sum of co- 
herent tunneling and sequential hoppings). Therefore, 
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dcphased transport is less sensitive to thermal fluctua- 
tions than coherent tunneling even if the 7r-orbital conju- 
gation is thermally disturbed. Another important feature 
in the log T deviation within dephased transport is that 
it decreases for longer molecules. This is simply caused 
by the fact that the number of paths through which the 
charge can be transmitted increases as the molecule gets 
longer. The averaging over the increased number of hop- 
ping pathways results in a smaller transmission disper- 
sion of dephased transport. 
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In summary, we have investigated electron transport 
in the presence of dephasing effects by analyzing both 
the coherent and incoherent process in molecular wires 
connected between two contacts. By employing the 
D'Amato-Pastawski model with first principle calcula- 
tions, two transport mechanisms could be examined un- 
der room temperature conditions where the structural 
conformations of the systems were explicitly treated. We 
were able to observe the tunncling-to-hopping transi- 
tion from exponential (coherent) to Ohmic (incoherent) 
regimes with molecular length using only one control pa- 
rameter. We also investigated the influence of thermal 
fluctuations and found that the transmission distribu- 
tions for tunneling transport are broader in comparison 
to those for the hopping mechanism. This highlights that 
the coherent component is rather more sensitive to struc- 
tural disorder than the incoherent one. Our simple ap- 
proach will provide a simulation platform for the inter- 
pretation of transport measurements in nanostructures. 
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